Nonlinear field theory with topological solitons: Skyrme models 
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£Nj ' In this talk, we give new insight into one of the best-known nonlinear field theories, 

the Skyrme model. We present some exact relevant solutions coming from different 



new versions (gauged BPS baby as well as vector BPS Skyrme models) giving rise 
to topological solitons, and highlighting the BPS character of the theory. 

PACS numbers: 

I. INTRODUCTION 



In 1834, the Scottish naval engineer J. Scott Russell saw in the Edinburgh canal a solitary 
wave for the first time. Nowadays, they are known as solitons, and among their main features 



we can point out that they are localized solutions which conserve their form, even after 
collisions. One specific kind of these objects are the interesting topological solitons. They 
are given by the collective excitations of the relevant degrees of freedom of the system and 
their existence and stability are caused by the global topological structure of the base and 
field space. 

It is here where the Skyrme model appears as a nonlinear field theory, introduced to 



describe nuclei, supporting topological solitons |10Ml2j. The main feature of this model is 
that it considers the primary fields as pions whereas nucleons and nuclei are described by 
collective nonlinear excitations of the degrees of freedom, i.e., by topological solitons. One 
of the most clever ideas of Skyrme was to identify the topological charge which arises in 
the model with the baryon number ensuring the latter to be conserved. Although this is 
a theory in 3+1 dimensions, due to its nonlinear character and complexity, sometimes it is 
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important to study its little brother in 2+1 too, the Baby Skyrme model 8|, |9), since it is 
simpler and can give us some clues to deal with the original model. This is an advice that we 
will follow in the present talk where two of the three different models are 2+1 dimensional 
theories. 

The Standard Skyrme Model consists of two terms: a quadratic term in first derivatives 
(the nonlinear sigma-model term), and a quartic one (the so-called Skyrme term). Both of 
them are necessary for the stability of the solutions from the point of view of the Derrick 
scaling, however, we can generalize the model with two more terms asking to have terms at 
most quadratic in time derivatives (so we have a standard Hamiltonian formulation). Thus, 
we can add a potential and a sextic term, so the more general model is given by 



with 



■5?Skr — + =^4 + =S?6 + =^0j (1) 



^ 2 = -^Tr(L^), ^ 4 = -^Tr ([L M , L v f) 
^ = X\ 4 Bl J? = -/iV, 



and 



B» = -—TrU^tfdvUtfdpUtfdvU), 

24:7T 2 

where B^ is the topological current and the degrees of freedom are the U fields (U G 
SU(2)). This effective field theory is related to QCD in the large N (number of colours) limit, 
which does not favour any particular term. Another important comment on this Skyrme 
model is that although originally it was thought for Nuclear and Particle Physics, now it 
has applications in other branches, as for instance the planar Skyrmions in ferromagnetic 
materials: a collective excitation of spins from the homogeneous case. 

Just to conclude this brief introduction we will present a key concept in our work, the 
BPS Bound, named after Bogomol'nyi, Prasad and Sommerfield. It consists in a lower 
energy bound of the system depending on its topology and not on the field configuration. 
Solutions saturating the bound (BPS solutions) fulfill a set of first order equations called 
BPS equations. Thus, BPS solutions imply a simplification of the field equations going 
from second to first order. Furthermore, when they saturate the bound, they minimize the 
energy, so BPS states are stable. An important characteristic for the Skyrme model is that 
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the bound is related to topological quantities as the winding number: E oc n. Unfortunately, 
the Standard Skyrme model cannot have solutions saturating the bound. One case where 

n 

this is possible is the BPS Skyrme Model [5 ! ] where we just take into account the sextic as 
well as the potential term: 

■^bps = =^6 + =^o- (2) 

This introduces the idea of infinite massive (quenched) pions since we forget about the 
kinetic term J2f 2 - This restricted Skyrme model solves old problems of hadrons and nuclei, 
and is being actively used in their phenomenology {4, 6, 7] 

With all these concepts in mind we are going to introduce three different new versions 
of the Skyrme model. In section 2, we will present a BPS baby model after introducing the 
gauge group U(l) with the usual Maxwell term for the gauge field. On the other hand, in 
section 3 we study the Skyrme field coupled to the vector meson u, firstly in 2+1 dimensions 
to finally reach the 3+1 dimensional case. 



II. GAUGED BPS BABY SKYRME MODEL 



The first model we are going to deal with is the gauged version of the BPS Skyrme Model 
above jlj] but in 2+1 dimensions, the Gauged BPS baby Skyrme model []]. The Lagrangian 
density is given by 

X = -^(Dj x D ^(j)) 2 - fi 2 V(n ■ $) - ^F a V (3) 

The degrees of freedom of the system are described by a vector of scalar fields = (0i, 2 , 3 ) 
with unit length <fi 2 = 1 and where the first term is quartic in first derivatives and at the 
same time is the square of the topological current, so it plays the role of the J^4 and =Sf 6 in 
the 3+1 version. Moreover, the J??2 term is not necessary for stability. On the other hand, 
since we have included the coupling to the gauge field, as well as the typical Maxwell term 
the usual derivatives have been replaced by the covariant ones: 

D a (j> = d a (f) + A a n x <p. (4) 

Notice that if we ask for finite energy field configurations in the full baby model (the j£?2 
term included), <fr(x,t) has to approach its vacuum value independently of the direction so 
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the base space M 2 can be compactified to a two-sphere S 2 , and the field is a map S 2 — > S 2 
characterized by an integer winding number or topological degree 



deg[0] = / d 2 x<p ■ di<j> x d 2 (j) — k, k e Z. 

4:71 



(5) 



Although this is not necessary for the BPS baby model we will assume it too because of 
consistency. 

The next step, in order to simplify things, will be to choose n = (0, 0, 1) and the standard 
static ansatz: 



^ sin f(r) cos n 92 ^ 



0(r, ip) 



sin/(r) sinn ip 
\ cos/(r) J 
A = A r = 0, A v = na(r), 

so the electric and magnetic fields are 

Ei = 0, B = d ± A 2 - d 2 A ± = 



(6) 
(7) 

(8) 



whereas the winding number is just deg[0] = n. It will be also useful to introduce the new 
variable and field 



y = r 2 /2, /i=-(l-cos/). 



(9) 



In this model, we will work with the so called old potential: V Q = 1 — 3 = 2h, a 2+1 version 
of the standard Skyrme potential used in 3+1 dimensions. Then, the boundary conditions 
are 

h(0) = 1 /(0) = 71, a(0) = (10) 
at the origin, and if we have compacton solutions 



at the radius, or 



a y (y = yo) = 0, h(y ) = h y (y ) = 



lim h(y) = 0, lim a y (y) = 



(11) 



(12) 



for exponentially or powerlike localized solutions. It can be shown that the old baby potential 
corresponds to compacton solutions. Although it seems we have five conditions notice that 
we are including a new parameter, the compacton radius yo. 
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Before going to the solutions, we will ask ourselves if a BPS bound can be found for this 
model. Fortunately, the answer is positive. Thus, regarding the static energy functional and 
after a non-trivial manipulation we arrive at 

E > E X 2 J d 2 xqW = An\n\E \ 2 (W') S 2 = 2n\n\E X 2 \W(fo = -1)|, (13) 

where we will call W the superpotential and n the winding number (making clear it is a 
topological bound). This bound will be saturated if the fields obey the BPS equations given 
by 

Q = W\ (14) 
B = -g 2 X 2 W, (15) 

where 

Q = i-D l $xD 2 $^q + e ij A i d j {n-$). (16) 

Let us comment on the superpotential W present in this BPS bound. It appears in the 
process of deriving the bound that W is not a free function but is restricted to globally 
existing (on the unit interval) solutions of the following superpotential equation 

X 2 W ,2 + g 2 X 4 W 2 = 2^ 2 V. (17) 

We give it this name since it is very similar to the superpotential equation of fake supergravity 



13l | . In fact, we can find this kind of equation in self-gravitating domain walls, where the 
terms have different signs, or in extremal supersymmetric black holes, where they enter with 
the same sign. It is also necessary to remind that in order to have BPS solitons the solution 
of this superpotential equation has to exist globally. Since it depends on the potential V, 
we will not have a solution for all V. For instance, if V has two vacua there are no solitons. 
It is also important to remark that in this model all solutions of the static equations are 
BPS (we will see it with the numerical solutions). 

In figures [1] and [2] we can see the numerical solutions we have found for different values 
of the parameter g (fixing A = n = 1). In the case of the value of /i 2 we just show cases 
corresponding to fi 2 = 0.1 since the situation is quite similar for other values. However, 
this does not happen when considering g. In figure [1] we see the solution with g — 0.1 
(similar situation has been found for g lower than this.), where the derivative of the field 
a corresponding to the magnetic field is here localized at the center. On the other hand, 
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(a) Function h with its derivative (dashed line). (b) Function a and its derivative (dashed line). 



FIG. 1: Solutions for g = 0.1 and [i 2 = 0.1. 



1.0 



-0.5 



-1.0 



0.2 



-OA ~ - - 0.6 



0.8 




-2.0 



(a) Function h with its derivative (dashed line). (b) Function a and its derivative (dashed line). 



FIG. 2: Solutions for g = 2 and /i 2 = 0.1. 



looking now at figure |2] we see a different situation, with the magnetic field being almost 
constant until the compacton boundary and the field h approaching the zero value suddenly. 
In fact, the higher the value of g, the more pronounced is the approach to zero. Finally, 
using the numerical solutions for several values of g we can see that indeed all soliton 
solutions are BPS as shown in figure [31 where the solid line represents the calculated bound 
E B = 2irn\ 2 \W(h = 1)|. 

III. VECTOR SKYRME MODELS 

In the model above we have just partially answered the question of how solitons interact 
with the lightest fields. We say partially because besides the U(l) Maxwell field we can 
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FIG. 3: Energy of static solutions (open circles) compared to the BPS bound (continuous 

line). 



consider another one to couple to solitons: the omega vector mesons, which after pions, 
are the lightest ones and their (phenomenologically motivated) coupling is just oj^B^. As 
we have commented before, the 2+1 version of the system can be a good starting point to 
then move on to the 3+1 dimensional case, and this will be the process we are going to 
follow here. Finally, since we are going to compare this vector model with its BPS brother 
without additional mesons, we will have to be aware of possible changes in the typical mass 
hierarchy: the BPS model is based on pions with infinite mass while here we are regarding 
that heavier mesons have finite mass. 



A. Vector BPS baby Skyrme model 

As it was commented above, the difference of this model with respect to the first one is 
just that now we couple the solitons to a vector meson instead of a gauge field. Then, the 
2+1 dimensional version of the Skyrme term is replaced here by the vector meson coupling 

if = -/rV(<fe) - \(d,u u - d u u,) 2 + \m 2 ujI + A't^B", (18) 
where _B M is the topological current given by 

B' l = ~^$-(djxd p $), (19) 

and we will use the generalized old baby potentials 

V=( i +^V. <*>!■ (20) 
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As before, we can see the vector field as a map from S 2 to S 2 after compactification of 
the base space. Thus, we can use the stereographic projection 

1 , ., . |2 



l + M 2 

with the axially symmetric ansatz: 



u + u, — i(u — u), \u\ —1), (21) 



M = /(r)e in< ^ u = u = u(r) = 0, (22) 



where n is the winding number. 

The boundary conditions to have solitons will be 



f(r = 


o) = 


= oo, 


f(r 


= R) 


= 0, 


cj'(r = 


= o) 


= 0, 


u)(r 


= R) 


= 0, 


f'(r = 


R) 


= 0, 


u'{r 


= R) 


= 



(23) 



where R will be finite or infinite depending on the potential. 

Studying now solutions for the model we will focus on the most interesting case, the 
potential with a = 2. It is so interesting because we have an analytical solution as well as 
a BPS bound. Then, defining the new field g = and variable x = r 2 /2 = -^p^y (where 
A is the initial A' up to a constant), we can solve the equation of motion for g 



1 K\(lVT+ffi) 

s{v) = KJm TTWv • (24) 

with K\ the modified Bessel function of the second type and f3 = whereas the equation 
for the vector meson is just 



ijj„ 



^TI9- (25) 



Finally, the BPS energy is 

S=-j(7)>W),^(l-|j + -), (26) 

which is nonlinear on n. Therefore, states with high n are unstable with respect to those 
with lower n. 

Other results either for massive vector mesons or massless ones are compared with those 
for the BPS baby model in figure HI The most important difference appears for a = 1 since 
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BPS baby 

a=l a=2 



t 

C-compacton 

exponential power-like 



Massless vector 

a=l a=2 



> A 

C-compacton 

logarithmical 

Massive vector 

a=l a=2 




4 

exponential power-like 



C-compacton 



FIG. 4: Comparing solutions of the BPS baby Skyrme model with those of the baby vector 

one. 

in this case, we get a C-compacton for the vector model. This means that at the boundary 
of the compacton there is a discontinuity, so that a Dirac delta is needed as a source term, 
what implies that the topological charge is screened. However, it is a positive difference 
since, because of the hierarchy problem, we have to remember that in the vector model the 
omega mesons are massless or have a finite mass whereas the lighter pions in nature are 
thought to be infinitely massive. Then, if the situation were the same the baby BPS Skyrme 
model would not be useful from the phenomenological point of view. 

In the same way we could think about sending the mass of the omega mesons to infinity. 
With this we will recover the original mass hierarchy since if we assume that pions have 
infinite mass the same should stand for omegas. Then, our Lagrangian density (Fl8|) would 
become 

% = -t?V(<h) + \M 2 ul + A'c^fl", (27) 
so the field equation for u mesons results in 

u » = -fi2 B » ■ (28) 
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Plugging this back into the Lagrangian density we arrive at 

J2f = - \^Bl (29) 

Thus, in this limit of infinite massive vector mesons we actually get the BPS baby Skyrme 
model. 



B. Vector BPS Skyrme model 

The next and natural step is to use the knowledge we have from the baby model for the 
3+1 dimensional version. Now, the corresponding Lagrangian density is jsj: 

& = -n 2 V(U, U ] ) - ~(9 M w„ - d^f + l -M 2 ul + X'u^, (30) 

where is the topological current given by 

B» = —^Tr(e^'"tfd v Urfd p Urfd ff U), (31) 

247T 

and as before we will use the generalized Skyrme potentials, which in 3+1 dimensions are: 

V = (^^-J = (1 - cos 0". (32) 

Furthermore, similar to the situation we had in 2+1 dimensions, the corresponding funda- 
mental field U G SU(2) defines now a map from the three-sphere in the compactified base 
space to the field space SU(2), a map which is characterized by a winding number. Then, a 
suitable parametrization for this U field is 

U = = cos f + i sin £ n ■ f , (33) 

with 

ft = ; — 7={u + u, —i(u — u), 1 — |m| 2 ). (34) 

1 + \u\ 2 

And as before, in order to simplify things, we will assume the static ansatz: 

Q 

u = u = u(r), 0Ji = 0, £ = f (r), u = tan -e in< ^. (35) 
Therefore, the boundary conditions for these new variables £ and u will be 

£( r = 0) = 7T, w(r = i2o) = 0, 
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BPS Skyrme 
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FIG. 5: Comparing solutions of the BPS Skyrme model with those of the vector one. 



co r (r = 0) = 0, u(r = R ) = 0, (36) 

where Rq can be finite or infinite. 

Focusing now on the solutions of the massless and massive vector models, we can compare 
them with the BPS Skyrme model as we have done for the baby case (see figure Then, 
we can see that again C-compactons with the source term appear in the vector models for 
a = 3/2, but not for a = 1 where there are no solitons. This is an important issue since as we 
have commented above they are completely different models because of the mass hierarchy. 
However, in contraposition to what we found in the massive vector baby model, here we 
do not have a BPS solution for a = 2, so why did we call this model BPS? The answer is 
that fortunately there exists a potential allowing BPS solitons both for massless and massive 
vectors, namely: 

Vbps = -(£ -cose sin 2 , (37) 
and the BPS energy calculated is 

E = ^-u,(0), (38) 
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which for the massless case is analytically solvable getting 

E=—^=\J-nV 2 . (39) 

On the other hand, from analytical solutions for finite mass mesons we saw that in this case 
the energy behaves as E ~ n. So for both cases we have the same problem of unstable 
solitons for high n 6//5 we had in the baby model. 



IV. CONCLUSIONS 



In this talk we have presented three new different BPS versions of the Skyrme model. In 
all of them we studied how solitons interact with the lightest fields: gauge and vector meson 
fields. Specially in the vector meson model we have made use of the simpler baby version of 
the model to help us with its generalization to the 3+1 case. In the case of the gauge field, 
this generalization is a more difficult task which lies outside the scope of this work. 

Another important and interesting characteristic of all models is that BPS bounds have 
been discovered for them. In the case of the gauge model we have seen that a rich structure 
supports this topological bound with a superpotential equation similar to that in fake su- 
pergravity. In this case, BPS solitons exist for monotonically growing potentials and even 
more, all the solutions are BPS. Regarding the vector meson models, this BPS bound only 
exists for a specific potential: a generalized old baby potential in 2+1 dimensions and a 
similar to the generalized Skyrme potential for the 3+1 case. 

Finally, it is necessary to emphasize that for some potentials the vector BPS models 
differ qualitatively from the BPS ones. It is a really important difference for the case of the 
standard pion mass potential, since in the BPS models we are assuming pions with infinite 
mass whereas in the vector ones we are adding omega mesons with finite mass, so we are 
changing the mass hierarchy. Therefore, if we had obtained different results, the BPS model 
would not be suitable from the phenomenological point of view. 
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